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Abstract This chapter covers the development of feedback control of superconduct¬ 
ing qubits using projective measurement and a discrete set of conditional actions, 
here referred to as digital feedback. We begin with an overview of the applications 
of digital feedback in quantum computing. We then introduce an implementation of 
high-fidelity projective measurement of superconducting qubits. This development 
lays the ground for closed-loop control based on the binary measurement result. 
A first application of digital feedback control is fast and deterministic qubit reset, 
allowing the repeated initialization of a qubit more than an order of magnitude 
faster than its relaxation rate. A second application employs feedback in a multi¬ 
qubit setting to convert the generation of entanglement by parity measurement from 
probabilistic to deterministic, targeting an entangled state with the desired parity 
every time. 


1 Digital feedback control in quantum computing 

Moving from proof-of-principle demonstrations of quantum gates and al¬ 
gorithms to fully fledged quantum hardware requires closing the loop between 
qubit measurement and control. There are different categories of quantum 
feedback control, depending on the type of measurement and feedback law 
used. For clarity, we first offer a classification of quantum feedback, similarly 
to that used in classical feedback. Then, we focus on the particular class of 
discrete-time, digital feedback. 

1.1 Classification of quantum feedback 

A first distinction is between continuous-time and discrete-time feedback. In 
the first case, measurement and control are continuous in time and concur¬ 
rent. An example is the stabilization of a qubit state using continuous partial 
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measurement, as discussed in Refs. [4, 5, 6, 7, 8]. In discrete-time feedback, 
instead, the conditional control is applied only after a measurement has been 
performed and processed. Here, we focus exclusively on discrete-time imple¬ 
mentations. This class can be further divided into two categories, analog and 
digital. We speak of analog feedback when the measurement result assumes 
a continuum of values and the feedback law is a continuous function of the 
result. An example is the experiment in Ref. 9, where the feedback controller 
first integrates the signal produced by a weak measurement and then applies 
the resulting coherent operation on the qubit. If the measurement has a finite 
set of possible results, instead, the possible feedback actions are also finite. 
We refer to this as digital feedback. The simplest example is qubit reset (sec¬ 
tion 4.2), in which a strong projective measurement collapses the qubit into 
either the ground or excited state. Here, a tt rotation brings the qubit to 
ground. Another interesting example is digital feedback using ancilla-based 
partial measurement [10, 11]. In this case, the measurement output is dis¬ 
crete, showing that partial measurement is not necessarily associated with 
analog feedback. In many applications, digital feedback forces determinism 
into one of the most controversial aspects of quantum mechanics, namely the 
measurement, whose result is intrinsically probabilistic. Looking at the action 
of digital feedback as a black box, we expect to see a definite output qubit 
state for a given input. In an ideal feedback scheme, measurement results and 
the conditioned operations vary at every run of the protocol, but the overall 
process is deterministic and the output state is always the same. For example, 
one can project a two-qubit superposition to a specific Bell state by combining 
a parity measurement with digital feedback (section 5). 

1.2 Protocols using digital feedback 

Several quantum information processing (QIP) protocols call for digital feed¬ 
back. One of the requirements for a quantum computer is efficient qubit ini¬ 
tialization [12]. Often, the steady state of a qubit does not correspond to 
a pure computational state jO) or jl), bur rather to a mixture of the two. 
Therefore, active initialization methods have been used in many QIP archi¬ 
tectures. Examples are laser or microwave initialization [13, 14, 15, 16] and 
initialization by control of the qubit relaxation rate [17, 18]. An alternative 
method, recently used with NV centers in diamond [19] and superconducting 
qubits (section 2), relies on projective measurement to initialize the qubits 
into a pure state. However, measurement alone cannot produce the desired 
state with certainty, since the measurement result is probabilistic. Closing a 
feedback loop based on this measurement turns the unwanted outcomes into 
the desired state. A qubit register must not only be initialized in a pure state 
at the beginning of computation, but often also during the computation. For 
example, performing multiple rounds of error correction is facilitated by re¬ 
setting ancilla qubits to their ground state after each parity check [20]. When 
using a qubit as a detector (e.g. of charge [21] or photon parity [22]), a fast 
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reset can be used to increase the sampling rate without keeping track of past 
measurement outcomes. 

Similarly, in the multi-qubit setting, digital feedback is key to turning 
measurement-based protocols from probabilistic to deterministic. An example 
is the generation of entanglement by parity measurement [23]. A parity meas¬ 
urement projects an initial maximal superposition state into an entangled 
state with a well-defined parity, i.e., with either even or odd total number of 
qubit excitations (section 5). However, once again, the outcome of the par¬ 
ity measurement is random. When running the protocol open-loop multiple 
times, the average final state has no specific parity and is unentangled. Only 
by forcing a definite parity using feedback can one generate a target entangled 
state deterministically. 

A variation of closed-loop control, named feedforward, applies control on 
qubits different from those measured. Feedforward schemes have already found 
application in quantum communication, where the main objective is the se¬ 
cure transmission of quantum information over a distance. In quantum tele¬ 
portation, a measurement on the Bell basis of two qubits projects a third 
qubit, at any distance, into the state of the first, to within a single-qubit ro¬ 
tation [24]. The measurement result determines which qubit rotation, if any, 
must be applied to teleport the original state. An extension of teleportation 
is entanglement swapping [24]. This protocol transfers entanglement to two 
qubits which never interact, and forms the basis for quantum repeaters [25], 
aiming to distribute entanglement across larger distances than allowed by a 
lossy communication channel. Here, measurement and feedback are used in 
every step to first purify [26] and then deterministically transfer entangled 
pairs to progressively farther nodes. 

In quantum computing, feedforward operations are at the basis of the first 
schemes devised to protect a qubit state from errors. The simplest protocol 
is the bit-flip code [27], which encodes the quantum state of one qubit into 
a an entangled state of three, and uses measurement of two-qubit operators 
(syndromes) in combination with feedback to correct for Gx (bit-flip) errors. 
Of similar structure is the phase-flip code, which protects against (phase- 
flip) errors. To protect against errors on any axis, the minimum size of the 
encoding is five qubits. In quantum error correction, projective measurement 
is more than a tool to detect discrete errors that have already occurred. In 
fact, the measurement serves to discretize the set of possible erorrs. Measuring 
the error syndromes forces one and only one of these errors to happen. This 
greatly simplifies the feedback step, which is now restricted to a finite set of 
correcting actions. 

While few-qubit error correction schemes are capable of correcting any 
single error, they require currently inaccessible measurement and gate fidel¬ 
ities. A more realistic approach is offered by topologically protected circuits 
such as surface codes [28], where errors as high as 1% are tolerated at the 
expense of a larger number of physical qubits required [29]. One cycle in a 
surface code, aimed at maintaining a logical state encoded in a square lattice 
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of qubits, includes the projective measurements of 4-qubit operators as error 
syndromes. When an error is detected on a data qubit, the corrective, coher¬ 
ent feedback operation is replaced by a change of sign in the operators for the 
following syndrome measurements involving that qubit. In other words, errors 
are kept track of by the classical controller rather than fixed [30, 31]. Bey¬ 
ond protecting a state from external perturbations, performing fault-tolerant 
quantum computing will require robustness to gate errors. In surface codes, 
single- and two-qubit gates on logical qubits are also based on projective 
measurements and in some cases require digital feedback to apply conditional 
rotations [32]. 

In addition to the gate model [12], digital feedback is central to the 
paradigm of measurement-based quantum computing [33]. In this approach, 
also called one-way computation, the initial state is an entangled state of 
a large number of qubits. All logical operations are performed by projective 
measurements. To make computation deterministic, feedback selects the meas¬ 
urement bases at each computational step, conditional on the measurement 
results. 


1.3 Experimental realizations of digital feedback 

Digital feedback has been employed for entanglement swapping with trapped 
ions [34] and for the unconditional teleportation of photonic [35], ionic [36, 37], 
and atomic [38, 39] qubits. In linear optics, feedforward has been used to im¬ 
plement segments of one-way quantum computing [40, 41, 42, 43, 44, 45] and 
for photon multiplexing [46]. In the solid state, the first approach to feed¬ 
back, of the analog type, was used to stabilize Rabi oscillations of a supercon¬ 
ducting qubit [47]. Soon after, digital feedback with high-fidelity projective 
measurement was introduced in the solid state, also using superconducting 
circuits [2, 48]. Recently, digital feedback has been extended to multi-qubit 
protocols with superconducting qubits (section 5 and Ref. 49) and NV centers 
in diamond [50]. 


1.4 Concepts in digital feedback 

The basic ingredients for a digital feedback loop are: 1) projective qubit 
readout and 2) control conditional on the measurement result (see Fig. la 
for the simplest single-qubit loop). The main challenge for (1) is to obtain a 
high-fidelity readout which is also nondemolition, thus leaving the qubits in a 
state consistent with the measurement result. A mismatch between measure¬ 
ment result and post-measurement qubit state will trigger the wrong feedback 
action (Fig. lb). The requirement for (2) is to minimize the time, or latency, 
between measurement and conditional action. Various sources contribute to 
latency: the time for the signal to travel from the sample to the feedback con¬ 
troller, the time for the feedback controller to process the signal and discretize 
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it, and the delay to the execution of the conditional qubit gates. If a trans¬ 
ition between levels occurs in one of the measured qubits during this interval, 
for instance because of spontaneous relaxation, its state becomes inconsistent 
with the chosen feedback action, resulting in the wrong final state (Fig. Ic). In 
feedforward protocols, such as error correction or teleportation, the feedback 
action is applied to data qubits, which are different from the measured ancilla 
qubits. In this case, the loop must also be fast compared to the coherence 
times of data qubits. 
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Figure 1. Concept of a single-qubit digital feedback loop and possible 
errors, a, The measurement is digitized into either H or L for qubit declared in 
|0) or |1). A different unitary rotation is applied for each result. Errors occurring in 
case of qubit relaxation between measurement and action (b) or wrong measurement 
assignment (c). Top (bottom) row indicates the actual qubit state corresponding to 
result H (L). 


The simplest example of digital feedback is single-qubit reset. Here, the 
qubit is projected by measurement onto |0) or |1) and, depending on the 
targeted state, a tt pulse is applied conditional on the measurement result. 
In this example, we consider the effect of the errors in Fig. 1 b,c, modeling 
the qubit as a classical three-level system, where the third level includes the 
possibility of transitions out of the qubit subspace. This is relevant in the case 
of transmon qubits with a sizeable steady-state excitation [2, 48]. We indicate 
with pfj the probability of obtaining the measurement result M with initial 
state |z) and post-measurement state \j). With Fij we indicate the transition 
rates from |i) to |j), and with r-ph the time between the end of measurement 
and the end of the conditional operation. For perfect pulses, the combined 
errors for initial state cos(0) |0) -I- sin(6*) |1) are, to first order: 

-Per7° = Poo + POI + FoiTFb, 

= Pn + PlO + Pl 2 + (Ao + A 2 )T’Fb, 

and weighted combinations thereof for other 9. A simple way to improve feed¬ 
back fidelity is to perform two cycles back to back. While the dominant error 
for 0 = 0 remains unchanged, for 0 = tt it decreases to + P 12 + A 2 'rFb- 
The second cycle compensates errors arising from relaxation to |0) between 
measurement and pulse in the first cycle. However, it does not correct for 
excitation from |1) to |2). For this reason, adding more cycles does not signi¬ 
ficantly reduce the error, unless the population in |2) is brought back to the 
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qubit subspace. This can be done [2, 48] by a deterministic tt pulse returning 
the population from |2) to |1), or with a more complex feedback loop capable 
of resolving and manipulating all three states. 


1.5 Closing the loop in cQED 

Until recently, the coherence times of superconducting qubits bottlenecked 
both achievable readout fidelity and required feedback speed. The devel¬ 
opment of circuit quantum electrodynamics [51, 52] with 3D cavities (3D 
cQED) ]53] constitutes a watershed. The new order of magnitude in qubit 
coherence times (> 10 /rs), combined with Josephson parametric amplifica¬ 
tion ]54, 55], allows projective readout with fidelities ~ 99% and realizing 
feedback control with off-the-shelf electronics. In the following section, we 
detail our implementation of high-fidelity projective readout of a transmon 
qubit in 3D cQED. We then shift focus to the real-time signal processing by 
the feedback controller, and on the resulting feedback action. 



Cavity 

1 )-t>^ 

Qubit Ampiification 


Figure 2. Simplified schematic of a single-qubit feedback loop in cQED. 

Upon application of a measurement tone at cuc, the signal Uut obtained from pro¬ 
cessing of the cavity output, carrying information on the qubit state, is input to 
the feedback controller and compared to a preset threshold Vth- If Vout > Uh (or 
Vth < Uh), the conditional rotation 9 {9') is applied to the qubit. 


2 High-fidelity projective readout of transmon qubits 

2.1 Experimental setup 

Our system consists of an A1 3D cavity enclosing two superconducting trans¬ 
mon qubits, labeled Qa and Qb, with transition frequencies a;A(B)/27r = 
5.606(5.327) GHz, relaxation times Tia(b) = 23 (27) fjs. The fundamental 
mode of the cavity (TElOl) resonates at lOrj^T^ = 6.548 GHz (for qubits 
in ground state) with K/27r = 430 kHz linewidth, and couples with gi/27r ~ 
75 MHz to both qubits. The dispersive shifts ]52] Xa(b)/'^ = —3.7 (—2.6) MHz, 
both large compared to K/27r, place the system in the strong dispersive regime 
of cQED ]56]. 
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Qubit readout in cQED typically exploits dispersive interaction with the 
cavity. A readout pulse is applied at or near resonance with the cavity, and 
a coherent state builds up in the cavity with amplitude and phase encod¬ 
ing the multi-qubit state [52, 57]. We optimize readout of Qa by injecting a 
microwave pulse through the cavity at wrf = oJr — Xa, the average of the res¬ 
onance frequencies corresponding to qubits in |00) and |01), with left (right) 
index denoting the state of Qb (Qa) (Figs. 3a,d). This choice maximizes the 
phase difference between the pointer coherent states. Homodyne detection of 
the output signal, itself proportional to the intra-cavity state, is overwhelmed 
by the noise added by the semiconductor amplifier (HEMT), precluding high- 
fidelity single-shot readout (Fig. 3c). We introduce a Josephson parametric 
amplifier (JPA) [54] at the front end of the amplification chain to boost the 
readout signal by exploiting the power-dependent phase of reflection at the 
JPA (see Figs. 3a,b). Depending on the qubit state, the weak signal transmit¬ 
ted through the cavity is either added to or subtracted from a much stronger 
pump tone incident on the JPA, allowing single-shot discrimination between 
the two cases (Fig. 3c). 


2.2 Characterization of JPA-backed qubit readout and 
initialization 

The ability to better discern the qubit states with the JPA-backed readout is 
quantified by collecting statistics of single-shot measurements. The sequence 
used to benchmark the readout includes two measurement pulses, Ma and 
Mb, each 700 ns long, with a central integration window of 300 ns (Fig. 4a). 
Immediately before Mb, a tt pulse is applied to Qa in half of the cases, in¬ 
verting the population of ground and excited state (Fig. 4b). We observe a 
dominant peak for each prepared state, accompanied by a smaller one over¬ 
lapping with the main peak of the other case. We hypothesize that the main 
peak centered at positive voltage corresponds to state jOO), and that the smal¬ 
ler peaks are due to residual qubit excitations, mixing the two distributions. 
To test this hypothesis, we first digitize the result of Ma with a threshold 
voltage Vth) chosen to maximize the contrast between the cumulative histo¬ 
grams for the two prepared states (Fig. 4c), and assign the value H{L) to the 
shots falling above (below) the threshold. Then we only keep the results of 
Mb corresponding to Ma = H. Indeed, we observe that postselecting 91% 
of the shots reduces the overlaps from ~ 6 to 2% and from ~ 9 to 1% in 
the H and L regions, respectively (Fig. 4d). This supports the hypothesis of 
partial qubit excitation in the steady state, lifted by restricting to a subset of 
measurements where Ma declares the register to be in jOO). Further evidence 
is obtained by observing that moving the threshold substantially decreases 
the fraction of postselected measurements without significantly improving the 
contrast [^ -1-0.1 (0.2)% keeping 85 (13)% of the shots] (Fig. 5b). Postselec¬ 
tion is effective at suppressing the residual excitation of any one qubit, since 
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Figure 3. JPA-backed dispersive transmon readout, a, Simplified diagram 
of the experimental setup, showing the input path for the readout signal carrying 
the information on the qubit state (RF, green) and the stronger, degenerate tone 
(Pump, grey) biasing the JPA. Both microwave tones are combined at the JPA and 
their sum is reflected with a phase dependent on the total power (b), amplifying the 
small signal. An additional tone (Null) is used to cancel any pump leakage into the 
cavity. The JPA is operated at the low-signal gain of ~ 25 dB and 2 MHz bandwidth. 

c, Scatter plot in the I — Q plane for sets of 500 single-shot measurements. Light red 
and blue: readout signal obtained with an RF tone probing the cavity for qubits in 
|00) and |01), respectively. Dark red and blue: the Pump tone is added to the RF. 

d, Spectroscopy of the cavity fundamental mode for qubits in |00) and |01). The RF 
frequency is chosen halfway between the two resonance peaks, giving the maximum 
phase contrast (163°, see inset on the right). Figure taken from Ref. 1. 


the |01) and |10) distributions are both highly separated from |00), and the 
probability that both qubits are excited is only ^ 0.2% . 

The performance of JPA-backed readout and the effect of initialization by 
measurement are quantified by the optimum readout contrast, defined as the 
maximum difference between the cumulative probabilities for the two prepared 
states (Fig. 5a). Without initialization, the use of the JPA gives an optimum 
contrast of 84.9%, a significant improvement over the 26% obtained without 
the pump tone. Comparing the deviations from unity contrast without and 
with initialization, we can extract the parameters for the error model shown in 
Fig. 5c. The model [1] takes into account the residual steady-state excitation 
of both qubits, found to be ~ 4.7% each, and the error probabilities for the 
qubits prepared in the four basis states. Although the projection into |00) 
occurs with 99.8 ± 0.1% fidelity, this probability is reduced to 98.8% in the 
time T = 2.4 fxs between Ma and Mb, chosen to fully deplete the cavity of 
photons before the tt pulse preceding Mb- We note that r could be reduced 
by increasing k by at least a factor of two without compromising Tia by the 
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Figure 4. Ground-state initialization by measurement, a, Pulse sequence 
used to distinguish between the qubit states (Mb), upon conditioning on the res¬ 
ult of an initialization measurement Ma- The sequence is repeated every 250 /rs. b, 
Histograms of 500 000 shots of Mb, without (red) and with (blue) inverting the pop¬ 
ulation of Qa with a TT pulse, c, Histograms of Ma, with Vth indicating the threshold 
voltage used to digitize the result, d, Mb conditioned on Ma = H to initialize the 
system in the ground state, suppressing the residual steady-state excitation. The 
conditioning threshold, selecting 91% of the shots, matches the value for optimum 
discrimination of the state of Qa- Figure taken from Ref. 1. 


Purcell effect [58]. By correcting for partial equlibration during r, we calculate 
an actual readout fidelity of 98.1 ± 0.3%. The remaining infidelity is mainly 
attributed to qubit relaxation during the integration window. 

As a test for readout fidelity, we performed single-shot measurements of a 
Rabi oscillation sequence applied to Qa, with variable amplitude of a resonant 
32 ns Gaussian pulse preceding Mb, and using ground-state initialization as 
described above (Fig. 5d). The density of discrete dots reflects the probabil¬ 
ity of measuring H or L depending on the prepared state. By averaging over 
^ 10 000 shots, we recover the sinusoidal Rabi oscillations without (white) and 
with (black) ground-state initialization. As expected, the peak-to-peak amp¬ 
litudes (85.2 and 96.7%, respectively) equal the optimum readout contrasts 
in Fig. 5a, within statistical error. 

2.3 Repeated quantum nondemolition measurements 

In an ideal projective measurement, there is a one-to-one relation between the 
outcome and the post-measurement state. We perform repeated measurements 
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Rabi rotation angle 9 (rad) 


Figure 5. Analysis of readout fidelity, a, Cumulative histograms for Mb 
without and with conditioning on Ma = H, obtained from data in Figs. 4c,d. 
The optimum threshold maximizing the contrast between the two prepared states is 
the same in both cases. Deviations of the outcome from the intended prepared state 
are: 8.9% (1.3%) for the ground state, 6.2% (2.1%) for the excited state without 
(with) conditioning. Therefore, initialization by measurement and postselection in¬ 
creases the readout contrast from 84.9% to 96.6%. b, Readout contrast (purple) and 
postselected fraction (black) as a function of Vth- c, Schematics of the readout error 
model, including the qubit populations in the steady state and at r = 2.4 /rs after 
Ma- Only the arrows corresponding to readout errors are shown, d, Rabi oscilla¬ 
tions of Qa without (empty) and with (full dots) initialization by measurement and 
postselection. In each case, data are taken by first digitizing 10 000 single shots of 
Mb into H or L, then averaging the results. Error bars on the average values are 
estimated from a subset of 175 measurements per point. For each angle, 7 randomly- 
chosen single-shot outcomes are also plotted (black dots at 0 or 1). The visibility of 
the averaged signal increases upon conditioning Mb on Ma ~ H. Figure adapted 
from Ref. 1. 
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to assess the nondemolition nature of the readout, following Refs. 59, 60. The 
correlation between two consecutive measurements, Mb and Me, is found to 
be independent of the initial state over a large range of Rabi rotation angles 
9 (see Fig. 6a). A decrease in the probabilities occurs when the chance to 
obtain a certain outcome on is low (for instance to measure Mb = H 
for a state close to |01)) and comparable to readout errors or to the par¬ 
tial recovery arising between Mb and Me- We extend the readout model 
of Fig. 5c to include the correlations between each outcome on Mb and 
the post-measurement state. The deviation of the asymptotic levels from 
unity, Ph\h = 0.99 and Pb\l = 0.89, is largely due to recovery during r, 
as demonstrated in Fig. 6b. From the model, we extrapolate the correla¬ 
tions for two adjacent measurements, Ph\h{'T' = 0) = 0.996 ± 0.001 and 
= 0) = 0.985 ± 0.002, corresponding to the probabilities that pre- 
and post-measurement state coincide. In the latter case, mismatches between 
the two outcomes are mainly due to qubit relaxation during Me- Multiple 
measurement pulses, as well as a long pulse, do not have a significant effect 
on the qubit state, supporting the nondemolition character of the readout at 
the chosen power. 



for two consecutive measurements Mb and Me, separated by r = 2.4 fis. Following 
an initial measurement pulse Ma used for initialization into |00) by the method 
described, a Rabi pulse with variable amplitude rotates Qa by an angle 9 along the 
i-axis of the Bloch sphere, preparing a state with R|oi) = sin^(6'/2). Red (blue): 
probability to measure Me = H (L) conditioned on having obtained the same result 
in Mb, as a function of the initial excitation of Qa. Error bars are the standard error 
obtained from 40 repetitions of the experiment, each one having a minimum of 250 
postselected shots per point. Deviations from an ideal projective measurement are 
due to the finite readout fidelity, and to partial recovery after Mb ■ The latter effect 
is shown in b, where the conditional probabilities converge to the unconditioned 
values, Ph = 0.91 and Pb = 0.09 for r Ti, in agreement with Fig. 4, taking into 
account relaxation between the tt pulse and Me- Error bars are smaller than the 
dot size. Figure taken from Ref. 1. 
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Josephson parametric amplification has become a standard technique for 
the high-fidelity readout of qubits in cQED. Since this experiment and the par¬ 
allel work in Ref. 61, projective readout of transmon [49, 62, 63] and flux [64] 
qubits has been performed using different varieties of Josephson junction- 
based amplifiers. The technology for these amplifiers continuously evolves to 
meet the needs of quantum circuits of growing complexity. One approach to 
high-fidelity readout of multiple qubits is to increase the amplifier bandwidth 
to include several resonators, each coupled to a distinct qubit [11]. Recent im¬ 
plementations in this direction included Josephson junctions in a transmission 
line [65], in low-Q resonators [66, 67], or in a circuit realizing a superconduct¬ 
ing low-inductance undulatory galvanometer (SLUG) [68]. Another approach 
for multi-qubit readout uses dedicated, on-chip Josephson bifurcation ampli¬ 
fiers [69]. 


3 Digital feedback controllers 

The input to a feedback loop in cQED is the homodyne signal obtained by 
amplification and demodulation of the qubit-dependent cavity transmission 
or reflection, as shown above. The response of the feedback controller is one 
or more qubit microwave pulses, which are generated and sent to the device 
(Fig. 2). This loop has a significant spatial extension, as the qubits sit in the 
coldest stage of a dilution refrigerator, while the feedback controller is at room 
temperature. A round trip involves 5 — 10 m of cable, which translates to a 
propagation time of 25 — 50 ns without accounting for delays due to filters and 
other microwave components. This physical limitation, which would require 
fast cryogenic electronics to be overcome, is only a small fraction of the total 
latency. A major source of delay is the processing time in the controller, 
combined with the generation or triggering of the microwave pulses for the 
conditional qubit rotations. The details of this process depend on the type of 
controller. We describe the first implementations below. 

The first realization of a digital feedback controller used commercial com¬ 
ponents for data sampling, processing, and conditional operations [2]. The core 
of the controller is an ADwin-Gold, a processor with a set of analog inputs 
and configurable analog and digital outputs. The ADwin samples the readout 
signal once, at a set delay following a trigger from an arbitrary waveform gen¬ 
erator (Tektronix AWG5014). This delay is optimized to maximize readout 
fidelity. A routine determines the optimum threshold for digitizing the readout 
signal. This voltage is then used to assign or L to the measurement. For the 
reset function in section 4.2, the ADwin triggers another arbitrary waveform 
generator (Tektronix AWG520) to produce a tt pulse when the outcome is L. 
Pulse timings and signal delays in the feedback cycle are illustrated in Fig. 7. 
The total time between start of the measurement and end of the feedback 
pulse is « 2.6 /rs, mainly limited by the processing time of the ADwin. 
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Figure 7. Digital feedback loop with an ADwin controller, a, Schematic 
of the feedback loop, consisting of an ADwin, sampling the signal, and a Textronlx 
AWG520, conditionally generating a qubit tt pulse, b. Timings of the feedback loop. 
The measurement pulse, here 400 ns long, reaches the cavity at t = 0. The ADwin, 
triggered by an AWG5014, measures one channel of the output homodyne signal 
(red: qubit in |0), blue: |1)), delayed by ~ 200 ns due to a low-pass filter at its 
input side. After comparison of the measured voltage a.t t = 0.6 /rs to the reference 
threshold, the AWG520 is conditionally triggered at t = 2.54 /rs, resulting in a tt 
pulse reaching the cavity at 2.62 ^s. Figure adapted from Ref. 2. 


To shorten the loop time, our second generation of digital feedback used 
a complex programmable logic device (CPLD, Altera MAX V), acquiring the 
signal following a 8-bit ADC, in place of the ADwin. This home-assembled 
feedback controller offers two advantages over the first: a programmable integ¬ 
ration window and a response time of 0.11 /rs (Fig. 8), an order of magnitude 
faster than the ADwin. As the feedback response time is now comparable or 
faster than the typical cavity decay time, active depletion of the cavity [70] 
will be required to take full advantage of the CPLD speed and further shorten 
the feedback loop. 

Further developments in the feedback controller replaced the CPLD with 
a field-programmable-gate-array (FPGA) to increase the on-board memory 
and enable more complex signal processing. For example, the FPGA allows 
different weights for the measurement record and maximal correlation with 
the qubit evolution. A FPGA-based controller has also been employed for 
digital feedback at ETH Zurich [49]. Recent developments at TU Delft and at 
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Time (gs) 

Figure 8. Digital feedback loop with a CPLD-based controller, a, Schemat¬ 
ics of the feedback loop, with an ADC and a CPLD (or FPGA) board replacing the 
ADwin in Fig. 7. b, Timings of the feedback loop. The CPLD samples the signal at 
every clock cycle (10 ns) and then integrates it over a window set by a marker of an 
AWG5014. The internal delay of the CPLD breaks down into the analog-to-digital 
conversion (60 ns) and the processing to compare the integrated signal to a calib¬ 
rated threshold, determining the binary output (50 ns). These timings are multiples 
of the clock (reduced to 4 ns in a recent FPGA-based implementation [71]). The 
total delay in Ref. [21] is increased to 2 /rs to let the cavity return to the ground 
state before the conditional tt pulse. 


Yale [72] include the pulse generation on a FPGA board, eliminating the need 
of an additional AWG. For comparison, Fig. 9 shows the setup that would 
be required for the 3-qubit repetition code [27] using our first generation of 
feedback (a) and the most recent one based on FPGAs (b, c). 


4 Fast qubit reset based on digital feedback 

4.1 Passive qubit initialization to steady state 

Our first application of feedback is qubit initialization, also known as re¬ 
set [12]. The ideal reset for QIP is deterministic (as opposed to heralded or 
postselected, see previous section) and fast compared to qubit coherence times. 
Obviously, the passive method of waiting several times Ti does not meet the 
speed requirement. Moreover, it can suffer from residual steady-state qubit 
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Figure 9. Hardware comparison for feedback control in the bit-flip code. 

The bit-flip code requires a two-bit digital feedback, acting on three qubits. Scaling 
the system in Fig. 7 would take an AWG520 for each qubit (a). A recent imple¬ 
mentation [71] performs readout signal processing and pulse generation on FPGA 
boards, resulting in the compact controller shown in b, c. 


excitations [1, 47, 61, 73], whose cause in cQED remains an active research 
area. The drawbacks of passive initialization are evident for our qubit, whose 
ground-state population i^o) evolves from states po and pi as shown in Fig. 10. 
With po and pi we indicate our closest realization 99% fidelity) of the ideal 
pure states |0) and |1). P|q) at variable time after preparation is obtained by 
comparing the average readout homodyne voltage to calibrated levels , as 
in standard three-level tomography [74, 75]. These populations dynamics are 
captured by a master equation model for a three-level system: 

(Ao 0 \ /^’|0)\ 

-f’li) ^ ^01 -Tw - A 2 Ai ^’ii) • (2) 

\P\2)) V 0 A2 -Ai/ V^|2>y 

The best fit to the data gives the qubit relaxation time A = 1/Ao = 
50 ± 2 ps and the asymptotic 15.5% residual total excitation. 
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Figure 10. Transmon equilibration to steady state. Time evolution of the 
ground-state population P|o) starting from states po and pi (notation defined in the 
text). Solid curves are the best fit to Eq. (2), giving the inverse transition rates 
= 50 ± 2 ps, T' 2 ’)^ = 20 ± 2 ps, = 324 ± 32 ps, = 111 ± 25 ps. From the 
steady-state solution, we extract residual excitations E|i)^as = 13.1 ± 0.8%, E| 2 >,as = 
2.4 ± 0.4%. Inset: steady-state population distribution (bars). Markers correspond 
to a Boltzmann distribution with best-fit temperature 127 mK, significantly higher 
than the dilution refrigerator base temperature (15 mK). Figure taken from Ref. 2. 


4.2 Qubit reset based on digital feedback 

Previous approaches to accelerate qubit equilibration include coupling to dis¬ 
sipative resonators [17] or two-level systems [18]. However, these are also 
susceptible to spurious excitation, potentially inhibiting complete qubit re¬ 
laxation. Feedback-based reset circumvents the equilibration problem by not 
relying on coupling to a dissipative medium. Rather, it works by projecting the 
qubit with a measurement {Mb, performed by the controller) and condition¬ 
ally applying a tt pulse to drive the qubit to a targeted basis state (Fig. 11). A 
final measurement {Me) determines the qubit state immediately afterwards. 
In both measurements, the result is digitized into levels H or L, associated 
with jo) and |1), respectively. The digitization threshold voltage Fth maxim¬ 
izes the readout fidelity at 99%. The tt pulse is conditioned on Mb = L to 
target |0) (scheme Fbo) or on Mb = H to target |1) (Fbi). In a QIP context, 
reset is typically used to reinitialize a qubit following measurement, when it 
is in a computational basis state. Therefore, to benchmark the reset protocol, 
we first quantify its action on pg and pi. This step is accomplished with a 
preliminary measurement Ma (initializing the qubit in pg by postselection), 
followed by a calibrated pulse resonant on the transmon 0 O 1 transition 
to prepare pi. The overlap of the Me histograms with the targeted region 
{H for Fbg and L for Fbi) averages at 96%, indicating the success of reset. 
Imperfections are more evident for 6 = tt and mainly due to equilibration of 
the transmon during the feedback loop. A detailed error analysis is presen- 
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-o- 9 = 0, uncond, 9 = 0, cond 

-1^ 9 = n, uncond, 9 = n, cond 



Homodyne voltage (V) 

Figure 11. Reset by measurement and feedback, a, Before feedback: histo¬ 
grams of 300 000 shots of Mb, with (squares) and without (circles) inverting the 
qubit population with a tt pulse. Each shot is obtained by averaging the homodyne 
voltage over the second half (200 ns) of a readout pulse. H and L denote the two 
possible outcomes of Mb, digitized with the threshold Vth, maximizing the contrast, 
analogously to section 2. Full (empty) dots indicate (no) postselection on Ma > Vps- 
This protocol is used to prepare po and pi, which are the input states for the feed¬ 
back sequences in b and c. b, After feedback: histograms of Me after applying the 
feedback protocol Fbo, which triggers a tt pulse when Mb = L. Using this feedback, 
~ 99% (92%) of measurements digitize to H for 0 = 0 (vr), respectively, c, Feedback 
with opposite logic Fbi preparing the excited state. In this case, ~ 98% (94%) of 
measurements digitize to L for 0 = 0 (tt). Figure taken from Ref. 2. 
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Figure 12. Deterministic reset from any qubit state. Ground-state popu¬ 
lation P|o) as a function of the initial state pe, prepared by coherent rotation after 
initialization in po, a-s in Fig. 11. The cases shown are: no feedback (circles), Fbo 
(squares), Fbi (diamonds), twice Fbo (upward triangles), and Fbo followed by Fbi 
(downward triangles). The vertical axis is calibrated with the average measurement 
outcome for the reference states po,pi, and corrected for imperfect state preparation. 
The curve with no feedback has a visibility of 99%, equal to the average preparation 
fidelity. Each experiment is averaged over 300 000 repetitions. Inset: error probabil¬ 
ities for two rounds of feedback, defined as where \t} £ {0,1} is the target 

state. The systematic ~ 0.3% difference between the two cases is attributed to er¬ 
ror in the tt pulse preceding the measurement of following Fbi. Curves: model 
including readout errors and equilibration (section 1.4). 


ted below. We emphasize that qubit initialization by postselection is here only 
used to prepare nearly pure states useful for characterizing the feedback-based 
reset, which is deterministic. 

4.3 Characterization of the reset protocol 

An ideal reset function prepares the same pure qubit state regardless of its 
input. To fully quantify the performance of our reset scheme, we measure its 
effect on our closest approximation to superposition states |0) = cos(0/2) |0)-|- 
sin(0/2) |1). Without feedback, P|o) is trivially a sinusoidal function of 9, with 
near unit contrast. Feedback highly suppresses the Rabi oscillation, with P|q) 
approaching the ideal value 1 (0) for Fbo (Fbi) for any input state. However, 
a dependence on 9 remains, with Perr = 1 — F’lo) for Fbo (1 ~ -P|i> for Fbi) 
ranging from 1.2% (1.4%) for 0 = 0 to 7.8% (8.4%) for 9 = ir. The remaining 
errors are discussed in section 1.4. From Eqs. (1), using the best-fit P^- and 
Tpb = 2.4 ns, errors due to equilibration sum to 0.7% (6.9%) for 0 = 0 (tt), 
while readout errors account for the remaining 0.4% (1.4%). In agreement 
with these values, concatenating two feedback cycles suppresses the error for 
0 = TT to 3.4%, while there is no benefit for 0 = 0 (1.3%). 
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Figure 13. Fast qubit reset. Initialization errors as a function of initialization 
time Tinit under looped execution of a simple experiment leaving the qubit ideally 
in |1) (a, measurement and tt pulse) or |0) (b, measurement only). Empty symbols: 
initialization by waiting (no feedback). Solid symbols: initialization by feedback, 
with three rounds of Fbo and a tt pulse on the 1 <->■ 2 transition. Two data sets 
correspond to two different cooldowns: the one corresponding to Figs. 10-12 (black) 
and a following one with improved thermalization (blue). Curves correspond to a 
master equation simulation assuming perfect pulses and measured transition rates 
Fij (dashed, no feedback; solid, triple Fbo with a tt pulse on 1 2). Feedback 

reset successfully bounds the otherwise exponential accruement of Perr in case a as 
Tinit —t 0. The reduction of Perr in b reflects the cooling of the transmon by feedback 
(see text for details). Figure adapted from Ref. 2. 


4.4 Speed-up enabled by fast reset 

The key advantage of reset by feedback is the ability to ready a qubit 
for further computation fast compared to coherence times available in 3D 
cQED [53, 76]. This will be important, for example, when refreshing ancilla 
qubits in multi-round error correction [20]. We now show that reset suppresses 
the accumulation of initialization error when a simple experiment is repeated 
with decreasing in-between time Tinit. The simple sequence in Fig. 13 emulates 
an algorithm that leaves the qubit in |1) [case (a)[ or |0) [case (b)[. A meas¬ 
urement pulse follows Tinit to quantify the initialization error Perr- Without 
feedback, Perr in case (a) grows exponentially as Tinit —0. This accruement of 
error, due to the rapid succession of tt pulses, would occur even at zero tem¬ 
perature, where residual excitation would vanish (i.e., Pi+i,i = 0), in which 
case Perr —> 50% as Tinit 0. In case (b), Perr matches the total steady-state 
excitation for all Tinit- Using feedback significantly improves initialization for 
both long and short Tinit- For Tinit 3> Tf, feedback suppresses Perr from the 
16% residual excitation to 3% (black symbols and curves)^, cooling the trans- 

^ We note that P|l> « P|2> = 1.6% is a non-thermal distribution. 
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mon. Crucially, unlike passive initialization, reset by feedback is also effective 
at short Tinit, where it limits the otherwise exponential accruement of error 
in (a), bounding Perr to an average of 3.5% over the two cases. Our scheme 
combines three rounds of Fbo with a pulse on the 10 2 transition before 
the final Fbo to partially counter leakage to the second excited state, which 
is the dominant error source [see Eq. (1)]. The remaining leakage is propor¬ 
tional to the average P|i), which slightly increases in a and decreases in b as 
Rnit —> 0. In a following cooldown, with improved thermalization and a faster 
feedback loop (Fig. 8), reset constrained Perr ^ 1% (blue), quoted as the fault- 
tolerance threshold for initialization in modern error correction schemes [29]. 
In addition to the near simultaneous implementation at ENS [48[, similar im¬ 
plementations of qubit reset have followed at Yale [72[ and at Raytheon BBN 
Technologies using a FPGA-based feedback controller. 

5 Deterministic entanglement by parity measurement 
and feedback 

In this section, we extend the use of digital feedback to a multi-qubit experi¬ 
ment, targeting the deterministic generation of entanglement by measurement. 
We first turn the cavity into a parity meter to measure the joint state of two 
coupled qubits. By carefully engineering the cavity-qubit dispersive shifts, we 
make the cavity transmission only sensitive to the excitation parity, but unable 
to distinguish states within each parity. Binning the final states on the parity 
result generates an entangled state in either case, with up to 88% fidelity to 
the closest Bell state. Integrating the demonstrated feedback control in the 
parity measurement, we turn the entanglement generation from probabilistic 
to deterministic. 

5.1 Two-qubit parity measurement 

In a two-qubit system, the ideal parity measurement transforms an un¬ 
entangled superposition state [■i/'®) = (|00) -I- jOl) -I- jlO) -I- |ll))/2 into Bell 
states 

\<P+) = ^(jOl) + jlO)) and = ^(jOO) + jll)) (3) 

for odd and even outcome, respectively. Beyond generating entanglement 
between non-interacting qubits [23, 77, 78, 79, 80[, parity measurements allow 
deterministic two-qubit gates [81, 82[ and play a key role as syndrome de¬ 
tectors in quantum error correction [24, 83[. A heralded parity measurement 
has been recently realized for nuclear spins in diamond [84[. By minimizing 
measurement-induced decoherence at the expense of single-shot fidelity, highly 
entangled states were generated with 3% success probability. Here, we real¬ 
ize the first solid-state parity meter that produces entanglement with unity 
probability. 
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Figure 14. Cavity-based two-qubit parity readout in cQED. a, Simplified 
diagram of the experimental setup. Single- and double-junction transmon qubits (Qa 
and Qb, respectively) dispersively couple to the fundamental mode of a 3D copper 
cavity enclosing them. Parity measurement is performed by homodyne detection of 
the qubit state-dependent cavity response [51] using a JPA [54]. Following further 
amplification at 4 K (HEMT) and room temperature, the signal is demodulated and 
integrated. A FPGA controller closes the feedback loop that achieves deterministic 
entanglement by parity measurement (Fig. 17). b, Matching of the dispersive cav¬ 
ity shifts realizing a parity measurement, c, Ensemble-averaged homodyne response 
(Vp) for qubits prepared in the four computational basis states, d, Curves: corres¬ 
ponding ensemble averages of the running integral (Vint) of (1V>) between fi = 0 
and = t. Single-shot histograms (5 000 counts each) of Vjnt are shown in 200 ns 
increments. Figure adapted from Ref. 2. 


5.2 Engineering the cavity as a parity meter 

Our parity meter realization exploits the dispersive regime [51] in two-qubit 
cQED. Qubit-state dependent shifts of a cavity resonance (here, the funda¬ 
mental of a 3D cavity enclosing transmon qubits Qa and Qb) allow joint qubit 
readout by homodyne detection of an applied microwave pulse transmitted 
through the cavity (Fig. 14a). The temporal average Vint of the homodyne re¬ 
sponse Vp{t) over the time interval [tj, tf] constitutes the measurement needle, 
with expectation value 

(Vint) =Tr (Op), 

where p is the two-qubit density matrix and the observable O has the general 
form 

O = l3o+ Pact^ + /Sect® -b PBAcrfa^. 

The coefficients /3o, Pa, Pb, and /3 ba depend on the strength Cp, frequency /p 
and duration rp of the measurement pulse, the cavity linewidth k, and the 
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frequency shifts 2xa and 2xb of the fundamental mode when Qa and Qb 
are individually excited from |0) to |1). The necessary condition for realizing 
a parity meter is ^a = /3b = 0 (^o constitutes a trivial offset). A simple 
approach [85, 86 ], pursued here, is to set /p to the average of the resonance 
frequencies for the four computational basis states \ij) {i,j G { 0 , 1 }) and to 
match xa = Xb- We engineer this matching by targeting specific qubit trans¬ 
ition frequencies /a and /b below and above the fundamental mode during 
fabrication and using an external magnetic field to fine-tune /b in situ. We 
align XA to xb to within ~ 0.06/t = 27r x 90 kHz (Fig. 14b). The ensemble- 
average (Hp) confirms nearly identical high response for odd-parity computa¬ 
tional states | 01 ) and | 10 ), and nearly identical low response for the even-parity 
100) and |11) (Fig. 14c). The transients observed are consistent with the in¬ 
dependently measured k, xa a-nd xb values, and the 4 MHz bandwidth of 
the JPA at the front end of the output amplification chain. Single-shot histo¬ 
grams (Fig. 14d) demonstrate the increasing ability of Vint to discern states of 
different parity as t{ grows (keeping ti = 0 ), and its inability to discriminate 
between states of the same parity. The histogram separations at t{ = 400 ns 
give |/3 a|, |/3b| < 0.02 |/3ba|- 

5.3 Two-qubit evolution during parity measurement 

Moving beyond the description of the measurement needle, we now investigate 
the collapse of the two-qubit state during parity measurement. We prepare the 
qubits in the maximal superposition state |V'°) = 5 (|00) -I- |01) -I- |10) -I- |11)), 
apply a parity measurement pulse for rp, and perform tomography of the final 
two-qubit density matrix p with and without conditioning on Vint (Fig. 15a). 
We choose a weak parity measurement pulse exciting hss = 2.5 intra-cavity 
photons on average in the steady-state, at resonance. A delay of 3.5 /k = 350 ns 
is inserted to deplete the cavity of photons before performing tomography. The 
tomographic joint readout is also carried out at /p, but with 14 dB higher 
power, at which the cavity response is weakly nonlinear and sensitive to both 
single-qubit terms and two-qubit correlations (/3 a ~ /3b ~ Pba, as required 
for tomographic reconstruction [87]. 

The ideal continuous parity measurement gradually suppresses the uncon¬ 
ditioned density matrix elements pij^ki = connecting states with 

different parity (either i ^ k or j ^ 1), and leaves all other coherences (off- 
diagonal terms) and all populations (diagonal terms) unchanged. The exper¬ 
imental tomography reveals the expected suppression of coherence between 
states of different parity (Figs. 15b,c). The temporal evolution of jpnpol) 
with near full suppression by rp = 400 ns, is quantitatively matched by a 
master-equation simulation of the two-qubit system. Tomography also un¬ 
veils a non-ideality: albeit more gradually, our parity measurement partially 
suppresses the absolute coherence between equal-parity states, jpoi.iol and 
|poo,ii|- The effect is also quantitatively captured by the model. Although 
intrinsic qubit decoherence contributes, the dominant mechanism is the dif- 
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Figure 15. Unconditioned two-qubit evolution under continuous par¬ 
ity measurement, a, Pulse sequence including preparation of the qubits in the 
maximal superposition state parity measurement and tomography 

of the final two-qubit state p using joint readout, b, Absolute coherences |pii,io|, 
|poi,io|, |poo,ii| following a parity measurement with variable duration rp. Free para¬ 
meters of the model are the steady-state photon number on resonance has = 2.5±0.1, 
the difference (xa — XB)/7r = 235 ± 4 kHz, and the absolute coherence values at 
rp = 0 to account for few-percent pulse errors in state preparation and tomography 
pre-rotations. Note that the frequency mismatch differs from that in Fig. 14b due 
to its sensitivity to measurement power, c, d, Extracted density matrices for rp = 0 
(c) and rp = 400 ns (d), by which time coherence across the parity subspaces (grey) 
is almost fully suppressed, while coherence persists within the odd-parity (orange) 
and even-parity (green) subspaces. Error bars correspond to the standard deviation 
of 15 repetitions. Figure taken from Ref. 3. 


ferent AC-Stark phase shift induced by intra-cavity photons on basis states 
of the same parity [86, 88, 89]. This phase shift has both deterministic and 
stochastic components, and the latter suppresses absolute coherence under 
ensemble averaging. We emphasize that this imperfection is technical rather 
than fundamental. It can be mitigated in the odd subspace by perfecting the 
matching of xb to XA, and in the even subspace by increasing xa,b/k 1.3 
in this experiment). 























24 


Diego Riste and Leonardo DiCarlo 



Figure 16. Probabilistic entanglement generation by postselected parity 
measurement, a, Histograms of Vint (rp = 300 ns) for the four computational 
states. The results are digitized into Mp = 1(—1) for Vp below (above) a chosen 
threshold, b, Parity readout fidelity Fp as a function of rp. We define Fp = 1 —Ce —to, 
with £e = p{Mp = — l|even) the readout error probability for a prepared even state, 
and similarly for £□. Data are corrected for residual qubit excitations (1 — 2 %). Error 
bars are smaller than the dot size. Model curves are obtained from 5 000 quantum 
trajectories for each initial state and rp, with quantum efficiencies rj — 0.25, 0.5, and 
1 for the readout amplification chain. No single value of rj matches the dependence 
of Fp on rp. We attribute this discrepancy to low-frequency fluctuations in the 
parametric amplifier bias point, not included in the model, c, Concurrence C of the 
two-qubit entangled state obtained by postselection on Mp = — 1 (orange) and on 
Mp — -1-1 (green squares). Empty symbols correspond to the threshold Vth that 
maximizes Fp, binning Psuccess ~ 50 % of the data into each case. Solid symbols 
correspond to a threshold Vth-(Vth+) for postselection on Mp = —1(-|-1), at which 
eo(ee) = 0.01. Concurrence is optimized at rp ~ 300 ns, where Paucceas ~ 20% in 
each case. We employ maximum-likelihood estimation [87] (MLE) to ensure physical 
density matrices, but concurrence values obtained with and without MLE differ by 
less than 3% over the full data set. d, e, State tomography conditioned on Vp > Vth- 
(d) and Vp < Vth+ (e), with rp = 300 ns, corresponding to the dark symbols in c. 
Figure taken from Ref. 3. 


5.4 Probabilistic entanglement by measurement and postselection 

The ability to discern parity subspaces while preserving coherence within 
each opens the door to generating entanglement by parity measurement on 
For every run of the sequence in Fig. 15, we discriminate using the 
threshold Vth that maximizes the parity measurement fidelity Fp (Fig. 16a). 
Assigning Mp = -1-1 (—1) to Vint below (above) Vth, we bisect the tomo¬ 
graphic measurements into two groups, and obtain the density matrix for 
each. We quantify the entanglement achieved in each case using concurrence 
C as the metric [90], which ranges from 0% for an unentangled state to 100% 
for a Bell state. As rp grows (Fig. 16b), the optimal balance between increas- 
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ing Fp at the cost of measurement-induced dephasing and intrinsic decoher¬ 
ence is reached at ~ 300 ns (Fig. 16c). Postselection on Mp = ±1 achieves 
C|Mp=-i = 45 ± 3% and C|Mp=-i-i = 17 ± 3%, with each case occurring with 
probability Psuccess ~ 50%. The higher performance for Mp = — 1 results from 
lower measurement-induced dephasing in the odd subspace, consistent with 
Fig. 15. 

The entanglement achieved by this probabilistic protocol can be increased 
with more stringent postselection. Setting a higher threshold Fth- achieves 
C|Mp=-i = 77 ± 2% but keeps Psuccess ~ 20% of runs. Analogously, us¬ 
ing Fth-i- achieves C|Mp=-i-i = 29 ± 4% with similar Psnccess (Figs. 16d, e). 
However, increasing C at the expense of reduced Psuccess is not evidently be¬ 
neficial for QIP. For the many tasks calling for maximally-entangled qubit 
pairs (ebits), one may use an optimized distillation protocol [90] to prepare 
one ebit from N = 1/Ej^f{p) pairs in a partially-entangled state p, where 
if At is the logarithmic negativity [90]. The efficiency S of ebit generation 
would he E = PsuccessEj\f (p). For postselection on Mp = —1, we calculate 
E = 0.31 ebits/run using Vth and E = 0.20 ebits/run using Vth-- Evidently, 
increasing entanglement at the expense of reducing Psuccess is counterproduct¬ 
ive in this context. 


5.5 Deterministic entanglement by measurement and feedback 

Motivated by the above observation, we finally demonstrate the use of feed¬ 
back control to transform entanglement by parity measurement from probab¬ 
ilistic to deterministic, i.e., Psuccess = 100%. While initial proposals in cQED 
focused on analog feedback schemes [91], here we adopt a digital strategy. Spe¬ 
cifically, we use our homebuilt programmable controller (section 3 to apply a 
TT pulse on Qa conditional on measuring Mp = -|-1 (using Vth, Fig. 17). In 
addition to switching the two-qubit parity, this pulse lets us choose which odd- 
parity Bell state to target by selecting the phase <p of the conditional pulse. To 
optimize deterministic entanglement, we need to maximize overlap to the same 
odd-parity Bell state for Mp = —1 (Fig. 17b) as for Mp = -fl (Fig. 17c). For 
the targeted state |^^), this requires cancelling the deterministic AC Stark 
phase pe = 0.737r accrued between jOO) and jll) when Mp = -|-1. This is 
accomplished by choosing p = {n — (pe)/2, which clearly maximizes the entan¬ 
glement obtained when no postselection on Mp is applied (Figs. 17c, d). The 
highest deterministic C = 34% achieved is lower than for our best probabilistic 
scheme, but the boost to Psuccess = 100% achieves a higher E = 0.41 ebits/run. 

A parallel development realized the probabilistic entanglement by meas¬ 
urement between two qubits in separate 3D cavities [92], establishing the first 
quantum connection between remote superconducting qubits. In another two- 
qubit, single-cavity system, feedback has been recently applied to enhance the 
fidelity of the generated entanglement [93]. Following the first realizations in 
3D cQED, parity measurements have been implemented using an ancillary 
qubit [62, 94] in 2D. Compared to the cavity-based scheme, the use of an 
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Figure 17. Deterministic entanglement generation using feedback, a, 

We close a digital feedback loop by triggering (via the FPGA) a tt pulse on 
Qa conditional on parity measurement result Mp = +1. This tt pulse switches 
the two-qubit parity from even to odd, and allows the deterministic targeting of 
1^^) = (|01) -I-110))/%/2. b, c, Parity measurement results Mp = —1 and Mp = -fl 
each occur with ~ 50% probability. The deterministic AC Stark phase acquired 
between |01) and |10) during parity measurement (due to residual mismatch between 
Xa and xb) is compensated by a global phase rotation in the tomography pulses. 
A different AC Stark phase is acquired between |00) and |11), resulting in the state 
shown in c, with the maximal overlap with even Bell state [|00)-|-exp(—ii/Pe) |ll)]/\/2 
at ifie = 0.7371. d, Generation rate of entanglement using feedback, as a function of 
the phase ip of the tt pulse. The deterministic entanglement generation efficiency 
outperforms the efficiencies obtained with postselection (Fig. 16). Error bars are 
the standard deviation of 7 repetitions of the experiment at each p. e, Full state 
tomography for deterministic entanglement [p = (tt — pe )/ 2 \, achieving fidelity 
(^"*"1 plfF*") = 66% to the targeted |^^), and concurrence C = 34%. Colored bars 
highlight the contribution from cases Mp = — 1 (orange) and Mp = -t-1 (green). 
Figure taken from Ref. 3. 


ancilla evades measurement-induced dephasing and is better suited to scaling 
to larger circuits. 


6 Conclusion 

We have presented the first implementation of digital feedback control in su¬ 
perconducting circuits, and its evolution to faster, simpler, and more config- 
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urable feedback loops. In particular, we showed the use of digital feedback for 
fast and deterministic qubit reset and for deterministic generation of entan¬ 
glement by parity measurement. Considering the vast range of applications 
for feedback in quantum computing, we hope that this development is just 
the start of an exciting new phase of measurement-assisted digital control in 
solid-state quantum information processing. 
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